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A new high-performance algorithm for calculating the spin- and parity-dependent shell model 
nuclear level densities using methods of statistical spectroscopy in the proton-neutron formalism was 
recently proposed. When used in valence spaces that cover more than one major harmonic oscillator 
shell, this algorithm mixes the genuine intrinsic states with spurious center-of-mass excitations. In 
this paper we present an advanced algorithm, based on the recently proposed statistical moments 
method, that eliminates the spurious states. Results for unnatural parity states of several sd-shell 
nuclei are presented and compared with those of exact shell model calculations and experimental 
data. 
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I. INTRODUCTION 

In this article we make a next step towards a reliable 
practical algorithm for calculating the level density in a 
finite many-body system with strong interaction between 
the constituents. Our primary object of applications is 
the atomic nucleus but the same techniques can be ap- 
plied to other mesoscopic systems, such as atoms in traps 
[l| and quantum dots 2]. Being a key element in the 
description of nuclear reactions and quantum transport 
in general, the many-body level densities are also inter- 
esting from the fundamental viewpoint since, at not very 
high excitation energy, the growth of the level density re- 
flects the interplay of interactions inside the system. The 
energy dependence of the level density may indicate the 
phase transformations smoothed out in finite systems: 
pairing quenching in nuclei [3, Q and small metallic par- 
ticles ;5], or magnetic effects in small quantum dots 
The correct reproduction of the level density serves as 
the first step to recognizing the regular or chaotic nature 
of spectral statistics in nuclei [1, 0, Q , complex atoms [9] 
and quantum dots ^IQj . 

The information on spin- and parity-dependent nuclear 
level densities (SPNLD) represents a critical ingredient 
for the theory of nuclear reactions, including those of 
astrophysical interest. For example, the routinely used 
Hauser-Feschbach approach [ll| requires the knowledge 
of nuclear level densities for certain quantum numbers 
J'^ of spin and parity in the window of excitation en- 
ergy around the particle threshold [l^, E^l • The nuclear 
technology requires the level density in the region of com- 
pound resonances close to the neutron separation energy. 
Therefore, a lot of effort has been invested in finding the 
accurate SPNLI), starting with the classical Fermi-gas 
approximations |14| - [la | and progressing to the sophisti- 
cated mean-field combinatorics [1, [3 and various 
shell model approaches with residual interactions in large 
valence spaces [l9l - [27j . 



Earlier we developed a strategy [T^ [26l - l29| of calcu- 
lating the shell model SPNLD using methods of statis- 
tical spectroscopy [sO, [HI . In the basic version of this 
approach, all basis states that can be built within a se- 
lected spherical single-particle valence space are taken 
into account. However, when the valence space spans 
more than one harmonic oscillator major shell, the shell 
model states include spurious excitations of the center- 
of-mass (CM). Apart from the level density problem, the 
correct separation of the spurious states is important for 
finding the physical response of the system to external 
fields, for example in the case of the excitation of isoscalar 
giant resonances. Over the years, shell model practition- 
ers developed techniques [33 that allow one to push these 
unwanted states to energies higher than those of interest 
for low-energy phenomena. These techniques involve (i) 
separation of basis states in the so-called Nhuj subspaces 
that can exactly factorize in the CM-excited and intrinsic 
states, and (ii) the addition to the nuclear Hamiltonian 
of a CM-part multiplied by a properly chosen positive 
constant. Unfortunately, the second ingredient leads to 
a multimodal distribution of levels and it is not appro- 
priate for the methods of statistical spectroscopy. 

In a recent letter [33j], we proposed a shell- model al- 
gorithm for removing the spurious CM contributions 
from the SPNLD that works if one knows the SPNLD 
for each Nhco subspace. These contributions are uni- 
modal and can be described using methods of statistical 
spectroscopy, provided that one can calculate the nec- 
essary moments of the Hamiltonian in Nhu) subspaces. 
In the present paper we formulate and utilize a high- 
performance algorithm that can calculate the configura- 
tion centroids and the widths of the Hamiltonian in Nhuj 
subspaces. This algorithm will be used for calculating the 
non-spurious level density for unnatural parity states of 
several srf-shcU nuclei. Comparison with the exact shell 
model results and/or experimental data will be also pre- 
sented. 
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FIG. 1: ^"Ne, Ihui, negative parity. Comparison of nuclear level densities from the exact shell model (stair lines) and from 
the moments method (straight lines). Dashed lines correspond to total densities with spurious states included; solid lines 
correspond to non-spurious densities without spurious states. Dotted lines present the spurious density. 



II. SPIN- AND PARITY-DEPENDENT 
CONFIGURATION MOMENTS METHOD 

We start closely following the approach proposed in 
Refs. [l^, H^, [131 . For reader's convenience, we repeat 
first the main equations of the moments method. Accord- 
ing to this approach, one can calculate the level density 
p as a function of excitation energy E in the following 
way: 

p{E,a)=Y,Do..-Ga,^iE), (1) 

K 

where a = {n, J,Tz,t:} is a set of quantum numbers, 
the total number of particles, n, total spin, J, isospin 
projection, T^, and parity, tt; for the level density, in 
contrast to the state density, the spin degeneracy 2J + 1 
is excluded. The sum over configurations k in Eq. ([1]) 
spans all possible (for the certain values of n,Tz, and tt) 
ways of distributing n particles over q spherical single- 
particle orbitals. Each configuration k is presented by a 
set of occupation numbers k — {ni, ri2, •.• , nq}, where rij 
is the number of particles occupying the spherical single- 
particle level j, J2j '^i = 

The energy dependence of the density p is expressed by 
finite-range Gaussian functions. Gam defined as in Ref. 

GaK{E) ^ G{E + Eg,s. — EaK,'^aK.), (2) 

s „ r exp {-x^/2a^) , \x\ < t] ■ a , . 
^{x,a)-C j Q , \x\>7j-a ' ^-^^ 

where the parameters Ea^ and ctqk will be defined later, 
Eg s. is the ground state energy, ij is the cut-off param- 
eter, and C is the normalization factor corresponding to 
the condition G{x, a)dx = 1. Although 77 can be 
treated as a free parameter, we know from the previous 
works (see for example [I^IH;!!^) that its optimal value 
is 77 ~ 3, which allows us to get a finite-range distribu- 
tion for the density and, practically, do not change each 



Gaussian contribution. Finally, the dimension in 
Eq. (HI) gives the correct normalization for each finite- 
range Gaussian being equal to the number of many-body 
states with given set of quantum numbers a that can be 
built for a given configuration k. 

The density distribution p{E), especially its low-energy 
part, is very sensitive to the choice of the ground state 
energy Egs . This origin of the energy scale is an external 
parameter for the moments method. To calculate it we 
need to use supplementary methods, such as the shell 
model. 

The parameters Ea,^ and ctqk in Eq. ^ are the fixed- 
J configuration centroids and widths. They essentially 
present the average energy and the standard deviation for 
the set of many-body states with a given set of quantum 
numbers a within a given configuration k. For a Hamil- 
tonian containing one- and (antisymmetrized) two-body 
parts of interaction, 

-ff = ^ Ciajfli + ^ X! Vijkiala^jaiak, (4) 

■i ijkl 

the fixed-J centroids and widths can be expressed in 
terms of the traces of the first and second power of this 
Hamiltonian [ii], Tr[H] and TiiH% fo r each configura- 
tion k: 

Eau = {H)aK, (5) 

aa. = Vm^. - {H)l,, (6) 

where 

(i/)„« = Tr("'=)[iJ]/i?„«, (7) 

{H^)a.^TT^"''^[H^]/Da.. (8) 

Every trace, such as Tr^"''-' [•••], contains the sum of all di- 
agonal matrix elements, '^(v, J\ - ■ ■ \ J), over all many- 
body states \i^,J) within given configuration k and with 
certain set of quantum numbers a. 
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As in our previous work [l9| we derive these traces for 
the fixed total spin projection J^, rather than for fixed 
J. Technicahy, it is much more easier to construct many- 
body states with a given total spin projection. The J- 
traces can be easily expressed through the J^-traces using 
the standard relation 

Tr(^^[---] = Tr('^^^[- -Tr('^^^[- (9) 

In Eq. (j9|) we omitted all quantum numbers, except for 
the projection and the total spin J. 

We will use the same label a to denote a set of quantum 
numbers that includes either the fixed Jz or the fixed J, 
keeping in mind that Eq. Q can always connect them. 
In every important case we will point out which set of 
quantum numbers was used. The expressions for the 
traces can be found in Ref. [s^], see also [l^ and ref- 
erences there. 



III. NEW FEATURES OF THE ALGORITHM 

To remove the center-of-mass spurious states from the 
level density we follow the approach suggested in [33| . 
This approach assumes the knowledge of the level den- 
sity within a restricted basis, and one needs to calculate 
partial densities p(E^ a, iV) for a given number N of ex- 
citations where the classification of states in terms 
of the harmonic oscillator Nhui excitations is applied. It 
has to be mentioned that such a splitting is clearly an 
approximation since it refers to the harmonic oscillator 
potential and the level scheme. The harmonic oscillator 
frequency being an auxiliary parameter has no strict 
a-priori definition. Nevertheless, as we will demonstrate, 
the suggested method results are independent of huj. 

According to this method a pure (without admixture 
of spurious states) level density p'"' {E, J, N) can be ex- 
pressed through the total (all states included) level den- 
sity p{E, J, N) for the same values of arguments and 
through the pure densities of lower N, Ref. fssj : 

pW(£;,j,iv) = p(i?,j,7V)- 

N Af.stcp 2 J+Jk 

-EE E P^"He,J\{N-K)). (10) 

K=l ,//f = J,„i„ J' = \J-Jk\ 

Here, for simplicity, we omitted all quantum numbers 
indicatting only total spin J representing the set a. In 
order to make these recursive equations complete we need 
a boundary condition that can be obtained from the Ohuj 
case which is free of spurious admixtures, 

p(o)(^,^,o) = />(£;, J, 0). (11) 

For example, if we are interested in the Ifiuj level density, 
we come to the following relation: 

J+i 

P^"He,J,1) = p{E,J,1)- PiE,J',0)- (12) 

J' = |,7-l| 



To calculate these partial level densities with the re- 
stricted values of excitation numbers N we need to 
slightly change the previous algorithm. By construction, 
each configuration k in Eq. ([T]) has a certain excitation 
number A^^ that can be defined as 

n 

Ar„ = ^i/a-A^o, (13) 

a=l 

where the sum runs over all n nucleons, Va is the excita- 
tion number of the single-particle level occupied by the 
nucleon a, and, for our convenience, we shift the result 
by the lowest value of the sum, A'o, so that N,^ starts 
with zero: N,^ — 0, 1, 2, • • •. Knowing A^^ we can restrict 
the sum in Eq. ([T]) by including only those configurations 
K which correspond to the desirable values of excitation 
numbers. 

For each configuration k we need to calculate the J- 
fixed width and centroid, cra„, Eau, defined by Eqs. ([5]) 
and There is no problem with the centroid calcu- 
lation, EaK ^ Tr^""-' [iJ] , since it remains unchanged, 
we just need to select the configurations k of our inter- 
est. The width calculation requires more attention. The 
problem is that the widths depend on the trace of the 
second power of the Hamiltonian, Tr^""-'[ff^], and if we 
want to restrict the basis we have to take care of the 
intermediate states — they also have to be restricted, 

Tr(-)[i72] = ^ J2 {X\H\p){p\H\X). (14) 

Here the sum over A includes all many-body states within 
the given configuration k and with certain set of quantum 
numbers a. In contrast to this, the sum over p includes 
all possible intermediate states regardless to configura- 
tions and quantum numbers, but only restricted by the 
desirable excitation numbers A^. In Eq. R repre- 

sents the set of all such many-body states defined by the 
allowed A^. In our previous algorithm, we could remove 
the intermediate sum using the completeness relation, 
J2fj. Im) (mI = 1; express the width in terms of the 
traces over single-particle excitations (I?— structures, as 
in Eq. (11) of Ref. [3). Now that is impossible, and we 
need to treat the situation difi'erently. 

One possible way to proceed here is still to follow the 
approach suggested in Ref. [13], where the "restricted" 
£)— structures were introduced, see Eqs. (11,12) in this 
reference. Here we propose an alternative solution, which 
may be easier to implement as an efficient computer al- 
gorithm. The single-particle part of the Hamiltonian (|4]) 
does not create any problems. Indeed the one-body oper- 
ators of the a|a.i type do not change the excitation num- 
ber A^, so that all states \p) in Eq. have the same 
Nfj^ as the states |A) all of which have the same excita- 
tion number defined by the configuration k, N\ = N,^. 
The two-body part of the Hamiltonian has the matrix 
elements 

- I Y.^^Mt^\aUhak\X) (15) 

ijkl 
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FIG. 2: ^^Na, Ifio;, negative parity. Comparison of nuclear level densities from the exact shell model (stair lines) and from 
the moments method (straight lines). Dashed lines correspond to total densities with spurious states included; solid lines 
correspond to non-spurious densities without spurious states. Right side of the figures presents the spurious density. 



which indeed can mix the states with different excita- 
tion numbers. Each member of the sum in Eq. (jlSp 
produces a certain change of the number of excitation 
quanta, AN^jki, 



AN, 



ijkl 



(16) 



where Vi is the excitation number of the single-particle 
level i. ft is important that ANijki reflects the internal 
property of the interaction operator ala^jaiak and does 
not depend on the many-body states. To deal with the 
sum over fi in Eq. (fH|) . we notice that A^^ = N\+ANijki, 
and therefore the restriction on N^^ can be reformu- 
lated in terms of the equivalent restriction on the single- 
particle levels, 



ti{N^,eR} ijkl 



E 

ijkl 



(17) 



where Ra means that the sum inside the parentheses is 
restricted by Nx + ANijki G R- 

In order to satisfy the condition of Eq. (17) we just 
need to keep those sets of the single-particle levels I, 
and k in (jlSp for which the excitation quanta of the |A) 
state plus the shift AiV, caused by the corresponding 
operator a|aja/afc, satisfy this restriction condition. For 
example, if we want to take into account only Ohuj and 
2hu) classes of states, and the |A) states belong to the, 
let's say, Ofioj class, then the allowed shifts of excitation 
quantum numbers will be AiV — 0, 2. All other terms 
in Eq. (jlSp must be ignored. Thus, we can use the old 
expressions for the ZJ-structures without changes (see Eq. 
(12) of Ref. [l^), which is an essential advantage, but 
it is necessary to enforce the appropriate restrictions on 
the sums over single-particle quantum numbers. 

The final result can be written as 



i i<j 

J2 2eiV,,.,D^^'^ 



-aK ' E ^m^qlqlD, 

{i<3mq<i) 



■Kj, q<l {4<j)/{q<0 



- E n.^^ 

i, (-3<0#i 



E2Vi - V, ( n[*J9l - 
^^liiq^lj]q \.^aK an 



I, (»<j)#9 



,(18) 



R2 



where i?i means that the sums inside the 
associated square bracket are restricted by 
(A''„ + Vi + Vj — Vq — vi) G R, and i?2 means that 
the sums in corresponding square bracket are restricted 
by (iV, + vq) e R. The trace D^l = Tr("")[a|a,] 
can be interpreted as a number of many-body states 
with fixed projection (if we consider Jz-traces) 
and the single-particle state i occupied, which can 
be constructed for the configuration k; the notations 
for more complex traces are ZJok 

D^a'^'^ = TT:^°"'\a\a]a\aqajai\, and so on. These D- 
structures were called propagation functions in Refs. 
p^ . [3^ . The detailed procedure for calculating these 
functions can be found in Ref. I19l. 
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In some applications, the restriction of the level den- 
sity to one class of excitations described by the excita- 
tion number iV, Eq. (10), might not be sufficient. For 
example, one could be interested in considering the class 
of (0 -I- 2)Uuj excitations. In those cases we select the 
maximum value iVmax of allowed excitations, and all the 
states, including intermediate ones, are restricted accord- 
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FIG. 3: ^^Mg, negative parity. Comparison of nuclear level densities between exact shell model (stair lines) and moments 
method (straight lines). Dashed lines correspond to total densities with spurious states. Solid lines correspond to non-spurious 
densities without spurious states. Right side of the figures present spurious density. 



ing to 



IV. RESULTS 



N < iV„ 



(19) 



Such a restriction allows us to fully take into account 
the interference between the states of different Nhw in 
Eq. (14). For example, if iVmax = 2, there are two 
possible classes of states, Ohuj and 2'huj, contributing to 
the width {Ihu! does not contribute because of oppo- 
site parity). The cross terms in Eq. (14), which are 
proportional to \{2huj\H\0hLL>)\'^ , are equally important 
along with the diagonal contributions \{2hLu\H\2huj)\'^ 
and \{Ofujj\H\Ohuj)\'^ . If one is interested in the density of 
unnatural parity states, one could consider only the Ihu; 
excitations. In these cases there are no cross terms, since 
the Ohuj states have opposite parity, thus one can directly 
use Eqs. (10) and (12). Finally, one should mention that 
the present algorithm was integrated in our hi ghly scal- 
able moments code that was described in Ref. |19| . 



As a first example we consider ^°Ne in the s — p — 
sd — p/-shell model space. For this space we use the 
Warburton-Brown (WBT) interaction [33 • All calcula- 
tions were done for Ihco subspace and for negative par- 
ity. Fig. 1 presents the comparison of the exact shell 
model level densities (stair-dashed and stair-solid lines) 
with the densities calculated using the moments methods 
(straight-dashed and straight-solid lines). The all dashed 
lines present the total densities including spurious states. 
The solid lines present the pure non-spurious densities. 
For the shell model calculations the spurious states were 
removed with the help of the Lawson method [32] that 
adds to the actual Hamiltonian a shifted center-of-mass 
Hamiltonian, Hcai, multiplied by a constant f3, 



H 



H' = H + P 



^^^^ - ^ 



(20) 



The additional term pushes the all spurious states up and 
leaves the non-spurious density at low-lying excitation 
energy. In our calculations we used /3 = 5. However, as 
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FIG. 4: ^*Si, all J, both parities. Comparison of nuclear level densities from experiment (stair lines) and from the moments 
method calculation (straight line). 



it was mentioned in the introduction, the Lawson recipe 
can not filter spurious states for the moments method. 
To get the straight-solid lines on Fig. 1 we use the re- 
cursive method introduced in Eqs. (10-12). Finally, the 
dotted lines present the spurious densities itself calcu- 
lated according to the second part of the right-hand side 
of Eqs. (10,12). To calculate the level density with the 
moments method we need to know the ground state en- 
ergy and the cut-off parameter rj. For ^"^Ne, the ground 
state energy £;g.s.(^°Ne) = -184.2 MeV was calculated 
with the help of the shell model, WBT interaction, OTiw 
subspace, and rj = 2.8. 

Figs. 2 and 3 present similar results for ^^Na, ^^Mg, 
in the s — p — sd — pf model space. The calculations also 
were done with the WBT interaction, for Ihuj subspace 
and negative parity. Stair lines refer to the shell model 
calculations, while straight lines present the results of the 
moments method. The only difference from Fig. 1 is the 
position of spurious states. The position of the spurious 
contribution to the total level density calculated with the 
shell model is naturally defined by the Lawson term and 
by parameter /3. As we can see from the Figures, (3 = 5 
is shifting all spurious states to the region of excitation 
energies of order of 140 MeV. The spurious states calcu- 
lated with the moments method are mostly located near 
the ground state energy. To compare the shape of the 
spurious part of the level density, we artificially shifted 
the results of the moments method by energy fHA MeV, 
which is 110 MeV for ^ = 22 and /3 = 5; after this, as 
seen from the Figures, the spurious states calculated with 
the shell model and Lawson term almost completely coin- 
cide with those calculated with the moments method and 
shifted afterwards. For all these cases we chose rj = 2.8 
and ground state energies Eg,s. = —201.2 MeV for ^^Na 
and Eg,s. = -202.3 MeV for '^'^Mg. 

Figures 4 and 5 present bigger cases of ^®A1 and ^®Si 
in the s — p — sd — pf model space, for both positive 
and negative parities. The dimensions are very large, 
and it is not practical to get the level densities with the 



shell model. Only the ground state energies can be cal- 
culated for these cases in OTioj subspace. With the WBT 
interaction, we got Eg,s. = -250.3 MeV for ^'^Al and 
Eg,s. = -285.0 MeV for ^*Si. Figures 4 and 5 show com- 
parisons of the level density calculated with the moments 
method, that is presented by straight lines, with experi- 
mental level densities presented by stair lines. There are 
two stair-like lines on each figure: the solid stair lines 
present an "optimistic" approach, when all experimen- 
tal levels with uncertain parity were counted; oppositely, 
the dashed stair lines present a "pessimistic" approach, 
when only experimental levels with defined parity were 
counted. 

It is needed to be mentioned that the real level densi- 
ties must be greater than those presented by the stair-like 
lines since it is possible that many levels were missed in 
experiment. In spite of the fact that the WBT interac- 
tion was not really tested in such big model spaces (it 
was specifically designed for A < 20), the agreement be- 
tween the calculations using the moments method and 
experimental data is quite remarkable. 



V. CONCLUSIONS AND OUTLOOK 

In conclusion, we developed a new high-performance 
algorithm to calculate the configuration centroids and 
widths of the nuclear Hamiltonian in Nhw subspaces 
built in a valence space. These first two moments can 
be used to calculate the SPNLD in the associated NTicj 
subspaces, which can be further mixed according to a 
recently proposed algorithm [s^] for extracting the non- 
spurious SPNLD. This strategy can be used to calculate 
the non-spurious shell-model level density for unnatural 
parity, a long-standing problem in nuclear structure. 

We tested our techniques by calculating the negative 
parity level density for several even-even sc?-nuclei, where 
the exact Ifiuj shell model diagonalization can be done 
and compared with the results of the newly proposed al- 
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FIG. 5: ^^Al, all J, both parities. Comparison of nuclear level densities from experiment (stair lines) and from the moments 
method calculation (straight line). 



gorithm. In all cases the SPNLD results of our moments 
method describe very well the results of the full shell 
model calculations. We also compared the results of our 
moments method with the available experimental data 
for ^^Al and ^®Si. For the states of positive parity, the 
effective interaction is well suited and our level density 
compares successfully with the experimental data. For 
the negative parity states, there is no well-adjusted Ihuj 
effective interaction for the middle of the sd-shell. Us- 
ing the WBT interaction that was tested up to about 
mass 20 [s^ we obtain a reasonable description of the 
experimental data. 

Certainly, there is a need of refined effective interac- 
tions for the unnatural parity states for the sd nuclei, and 
our techniques could help improving their quality. The 
algorithm suggested in this article can be used to reliably 
predict for the first time the SPNLD for a large number 
of unstable nuclei relevant for the rp-process. The next 



steps could be directed to the more complicated cases, 
where the Nhuj space is incomplete. This will allow us 
to reformulate the approach for the realistic mean-field 
potentials when it will be possible to exclude the stan- 
dard reference to harmonic oscillator and open the broad 
field of problems related to the cross sections and reaction 
rates. In parallel, the deep question should be addressed 
of the infiuence of continuum effects [s^ on the density of 
levels which are in fact resonance quasistationary states. 
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